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A linear system x = Ax, A G M"^", x e M", with rk^ = n~ 1, has a 
one-dimensional center manifold E'^ — {v G R" : Av — 0}- If a differential 
^0 I equation x — f{x) has a one-dimensional center manifold at an equi- 

■ librium x* then i?'^ is tangential to with A ~ Df(x*) and for stability 

of it is necessary that A has no spectrum in C+, i.e. if A is symmetric, 
it has to be negative semi-definite. 

We establish a graph theoretical approach to characterize semi-definiteness. 
Using spanning trees for the graph corresponding to A, we formulate meso- 
scale conditions with certain principal minors of A which are necessary for 
semi-definiteness. We illustrate these results by the example of the Ku- 
J> " ramoto model of coupled oscillators. 

^: 

m 

cn; 1 Introduction 

\ Center manifolds of differential equations x = f{x) play an important role, e.g. in 

bifurcation theory |6j |17j or, more specifically, if certain quantities or invariants are 
preserved, e.g. for mechanical systems with rigid-body modes, in chemical kinetics, 
compartmental modeling, and population dynamics, oftentimes these center mani- 
\ folds are one-dimensional and consist of a continuum of equilibria (see e.g. [3l |4j and 

the references therein). If a center manifold is one-dimensional and contains an 
equilibrium x* then ikDf{x*) = n ~ 1. For to be stable it is then necessary 
that A := Df{x*) has no eigenvalue with positive real part, i.e. A has to be nega- 
tive semi-definite, since we assume the Jacobian A to be symmetric. The assumption 
A = is e.g. satisfied if the coupling of the differential equations Xi = f{xi, . . . , a;„). 
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? = 1, . . . , n, is symmetric. W.l.o.g. one can assume that A = (aij) has row sum zero, 
i.e. '^ij — for i — 1, . . . , n. This is a consequence of the fact that the unit 

vectors vi e E" = {v G M" : Aw = 0} resp. wi = (f > ■ ■ ■ . f) £ can be 

extended to orthogonal bases {vi, . . . and {wi, . . . ,Wn} of M", i.e. the matrices 
V = {vi \ . . .\vn) and W = {wil . . . are orthogonaL The orthogonal transforma- 
tion y = WV^x preserves stability and transforms x — Ax and its center manifold E"^ 
into y = WV^ AVW^y with corresponding center manifold {awi £ M" : a G M}, i.e. 
Awi = 0. 

In the present paper we discuss semi-definiteness of a matrix A G R"^" by Sylvester's 
criterion. Our results provide a rigorous proof of a recently formulated stability condi- 
tion [7] and show that utilizing spanning trees of the Coates graph of A is a promising 
approach for identifying obstructions to the stability of one-dimensional center mani- 
folds, or formulated in a more application-oriented language, the impact of meso-scale 
structural properties on the dynamics of complex networks. As a motivational exam- 
ple consider the Kuramoto model of a network of coupled oscillators (see [TB] for an 
excellent survey) 

9, ^ujr+Y, B,, sin(0, -9{) , 'ii£{l...N} (1) 

on the A^-torus (5^)^, with = M/(27rZ) ~ [0, 27r) denoting the circle. Here, 

6i G and uji G denote the phase and the intrinsic frequency of node i, while 
(By ) G M^""^ is the symmetric weight matrix representing an undirected, weighted 
interaction network. Two oscillators «, j are thus connected if Bij = Bji ^ 0. Note 
that, if for a solution 6 = {9i, . . . , 6^) : E — > we sum up all equations in (IT|), then 
Eti S^{t) = Eti + Eti Ej^, B,j sm{9,{t) - e,{t)). Since B,, = it follows 
that Y.ti Y.J^^ B^J sm{ej{t) - e,{t)) = for t G K and hence e,{t) = Y.ti 

The Kuramoto model provides an example for studying synchronisation in networks 
and is widely accepted as a simple model of continous dynamics on networks. To see 
that the following three statements are equivalent 

(i) 9* is a phase-locked solution, i.e. 9*{t) = 9*{t) for all t G M, i, j G {1, . . . , N}, 

(ii) 9*{- + uj{T)) is a phase-locked solution for each io(T) G M, 

(iii) 9*{t) ~ fi for all t G M, G {1, . . . , N}, with the mean frequency D, := jj 

note that (iii) ^ (ii) => (i) is obvious, and (i) ^ (iii) follows from N9*{t) = 

^iLi^ii^) = J2iLi^i- ^ consequence, every phase- locked solution can be re- 
presented as t ^ 9*{t + uj{T)) = 9*{uj{T)) + {nt,...,nt) for an uj{T) G M and a 
phase-locked solution 9*. 

The Kuramoto model ([1]) can also be viewed in a 'rotating frame', more precisely, in 
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the new coordinates Xi{t) = 9i{t) — f2i. The transformed model 

Xi = LUi — n + Bij sm{xj ~ Xi) ^ Vi £ {1 . . . (2) 

has mean frequency equal to zero and every phase-locked solution 9* (• + w(T)), ijj{T) £ 
R, of (P is transformed to an equilibrium ^('^(^^ of © with x^"^*^'^ ^ e*{uj{T)), 
giving rise to a 1-dimensional manifold TV = {a;^"*^^-*^ £ : uj{T) £ M} consisting 
of a continuum of equilibria of ([2]). The corresponding manifold A4 — {9*{t) £ : 
i £ R and 0* is a phase- locked solution of ([T])}, consisting of all orbits of phase-locked 
solutions of ([TJ , is called stable with asymptotic phase if there exists a neighborhood 
U of A4 such that for any solution 6 starting in U, i.e. 6{0) £ U, there exists a phase- 
locked solution e* in 7W, i.e. 9*{t) £ X for t £ R, with limf^^o (0) = 0, 
where d{-,-) denotes the distance between two points on the torus . 

The linearization of ^ at an equilibrium £ A/" is given by i = Ax with 

Bijvcos(.^Ari) \ 

B2N C0s{(pN2) 

\ BniCOs{iPin) Bn2COs{iP2n) ••• ~J2j=£N ^Nj COs{ipjN)/ 

(3) 

where (pji :— 9*{lu{T)) - 9*{lu{T)) is independent of lo{T) £ R and the phase-locked 
solution 9*. Note that A has row sum zero and therefore is an eigenvalue. The 
following theorem characterizes the fact that the manifold A4 is stable with asymptotic 
phase by properties of ([2|) on and in the vicinity of TV. For an alternative approach 
using set-valued Lyapunov functions see e.g. [12] and the references therein. 

Theorem 1. Assume that 9* is a phase-locked solution of ([I} and rank A = N — 1. 
Then the following three statements are equivalent: 

(i) 9* is orbitally stable in the sense that M. is a stable manifold of phase-locked 
solutions of ([l} with asymptotic phase. 

(a) J\f is a stable center manifold of equilibria of ^ with asymptotic phase, i.e. 
there exists a neighborhood V of M such that solutions starting in V converge 
for t ^ oo to an equilibrium in J\f. 

(Hi) The manifold M of equilibria of ^ satisfies the following two conditions: 

(a) linear stability: A has no eigenvalues with positive real part, i.e. A is nega- 
tive semi-definite. 

(b) global attractivity: there exists a neighborhood V of Af, such that solutions 
.starting in V converge for t ^ oo to J\f. 

Proof. (Hi) (ii). This follows from 3, Proposition 4.1 & proof of Theorem 4.1]. 



A = 



B21 C0S(¥312) 



B12 COs{ip2l) 
-Ej^2 B2jCOs{(pj2) 
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{ii) => (iii). (iii)(b) follows directly from (ii). To prove (iii)(a), assume that A has 
an eigenvalue with positive real part. Then the linearization x = Ax has an unstable 
subspace giving rise to unstable manifolds for every equilibrium x'-"^"^-'-' of which 
contradicts the fact that every solution in V converges to an equilibrium in M for 
t ^ oo. 

{{) <^ {ii). The transformation Xi{t) = Oi{t) — fit maps M io N and preserves the 
property of asymptotic phase of corresponding solutions. □ 

Our main aim for the Kuramoto model in this paper is to 

(i) reveal the dependence of the stability of the continuum of equilibria on certain 
topological properties of the coupling network, 

(ii) efficiently evaluate algebraic stability criteria. 

To achieve this goal, we discuss a characterization of condition (iii) (a) of Theorem [1] 
which shows that semi-definiteness is a necessary condition for stability. 

Semi-definiteness can be characterized in terms of the determinants of the submatrices 
(minors) of a matrix via Sylvester's criterion. Stability analysis by means of Sylvester's 
criterion is well-known in control theory and has been applied to problems of different 
fields from fiuid- and thermodynamics to offshore engineering and social networks 
(see [8 and the references therein). From a modern perspective the application of 
Sylvester's criterion is very appealing as it can reveal stability criteria on different 
scales. The smallest minor considered by Sylvester's criterion is just a diagonal element 
of the Jacobian matrix and thus poses a condition on the internal dynamics on a single 
vertex. The subsequent minors pose conditions on subgraphs such as vertex pairs, 
triplets, and so on. Sylvester's criterion can thereby contribute to the search of meso- 
scale structural properties that have an impact on the global dynamics of networks, 
which is becoming a current topic in physics [2 I10| . However, due to the rapidly 
growing complexity of the minors to be evaluated, the application of the criterion 
was previously limited to systems with few degrees of freedom. We present a way to 
apply Sylvester's criterion to systems with many degrees of freedom. Generalizing the 
classical Kirchhoff Theorem, we obtain a characterization of the minors by properties 
of forests in the coupling graph of our system. 

We recall some basic notions from graph theory. A weighted graph G is a triple 
(y, E, ui) consisting of a finite non-empty set of vertices V, a set of edges E C l/^^) 
consisting of unordered pairs of vertices and a weight function cj : — > M \ {0}. 

A walk of length £ from ig £ V to i£ V is a sequence of vertices io, . . . ,ie such that 
{ik, ik+i} G E for k £ {1, — 1}. A cycle is a walk for which iq — ig and ik ^ iy 
for < /c < fc'. A graph is connected if for each pair of vertices (i,j) there is a walk 
from i to j. A forest is a graph without cycles, a tree is a connected forest. 



4 



A graph G' = {V',E',uj') is a subgraph of G = {V,E,uj) if V C V, E' C E and 
Lu' = LU\E'- If V' — V the subgraph is spanning. A spanning tree of a graph is a 
subgraph that is spanning and a tree. 

We will identify subsets E' C E with the induced subgraph G' — {V , E' ,uj') of G, 
where V ^ {i e V : 3j e V with {i,j} G E'} and cj' = cji^;/. 

A cut of the graph G given by a partition Vi , V2 of the vertex set is the set of edges 
E(Vi, V2) between Vi and V2. 

The weighted adjacency matrix A{G) G R"^" of a graph G {V,E) with F = 
{1, . . . , n} is defined by 



A{G\ 



\w{i,j) \i{i,j}&E 
otherwise 



For a symmetric matrix A — {aij) G M"^" the Coates-graph Q{A) — {V,E,oj) is 
defined as the weighted undirected graph with vertices V — edges E = 

aij 7^ 0} and weight function a;({i,j}) := aij. On the other hand, for each 
undirected graph its associated adjacency matrix is symmetric. 

The Laplace-matrix C{G) of a graph G is defined by C{G) = V(G)— A{G) where 2?(G) 
is the diagonal degree matrix 'D{G) :— {A{G)l)i. Take for example G with 



A{G) 



/O 

1 
2 


V-3 



-3\ 

-2 



k 

1 

1 1 
-2 1 / 



, then 



£(G) 




V 3 



3 \ 

-1 2 

2 -1 

-1 -4/ 



Remark 2. ('aj If loops in G are removed, then C{G) does not change. 

(h) If G has an adjacency matrix -4(G) which has row sum zero, then C{G) — —■4(G). 



The structure of the paper is as follows. In Section 2 we recall Sylvester's criterion 
in Theorem |21 Remark 2] emphasizes with an example the fact that testing for semi- 
definiteness is fundamentally more costly than checking definiteness. The main result 
of this section, Theorem[5l is a version of Sylvester's criterion for zero-row-sum matrices 
of maximal rank. Section 3 is devoted to the Matrix Tree Theorem [7] for principal 
minors. Its proof is split in three parts fLemmas [51 [51 and [TU|) . Section 4 contains 
Lemma [TT] on how to cut a forest, which is used to prove the combinatorial identity in 
Theorem 1121 Section 5 finally proves the existence of a positive spanning tree in the 
connected components of Coates graphs of negative semi-definite matrices and gives 
some applications of this theorem. 



5 



2 Definiteness and Sylvester's criterion 



In this section we recall Sylvester's criterion and adapt it to the case of positive semi- 
definite matrices with zero row sum and maximal rank. 

Let p,qeN,S C {!,..., p},T C {!,..., q} with IS*] = |T| and B e We denote 

the submatrix of B defined by S and T by Bs.t- The determinant of this submatrix 
is called the minor corresponding to S and T for which we write [_B]5_t. 

If i? is a square matrix, we call i?s,s a principal submatrix and its determinant a 
principal minor. If additionally S — {1, . . . , fc} for some fc £ N, we call it a leading 
principal submatrix /minor. 

RecaU that i e R"^" is called 

• positive (semi) definite iff G M" : Lv > [v^ Lv > 0), 

• negative (semi) definite iff Vw G M" : Lv < {v'^ Lv < 0). 

Theorem 3 (Sylvester's criterion, '181, Theorem 7.2]). ^ symmetric matrix L g M"^" 
is positive definite iff every leading principal minor of L is positive. L is positive 
semi- definite iff every principal minor of L is non-negative, i.e. [L]s,s ^ for all 
SC{l,...,n}. 

Remark 4 (Checking for semi-definiteness). Note that for definiteness one has to 
check the .sign of only n minors of L, while one has to check in general the sign of all 
2" principal minors to test for semi-definiteness. Even if we know that is a simple 
eigenvalue of L, Q{L) is connected and L has eigenvector 1, it is not enough to check 
the leading principal minors for non- negativity to ensure semi-definiteness. Consider 
the following matrix: 




The corresponding graph is connected, the corresponding characteristic polynomial is 
+ 2x^ — Ax'^ — Ax (hence zero is a simple eigenvalue), the leading principal minors 
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[C]{1,2,3},{1,2,3} 



1 

-1 1 

1 1 -2 



= 1 




-1 







[C]{1,2},{1,2} 

[C]{iUi} = |0|=0- 
But C is not positive semi-definite since for v = (1, 1, 0, 0)^ we have v'^Cv 



-1 < 0. 



The following theorem is a version of Sylvester's criterion for zero-row-sum matrices. 

Theorem 5 (Sylvester's criterion for zero-row-sum matrices of maximal rank). Let 
L G M"^" be a symmetric matrix with zero row sum. Then the following are equivalent: 

(i) L is positive semi-definite and has rank n — 1 

(i^) [L\s,s>() for alls ^{\,...,n} 

(Hi) Ls^s "is positive definite for all S C {1, . . . , n} 

M [L]s,s > for all S ^ {I, . . . ,k}, k e {1, . . . ,n ~ 1} 

(v) Ls^s is positive definite for S = {1, . . . , n — 1} 



Proof. The equivalences (ii) (Hi) and (iv) <=f- (v) follow directly from Sylvester's 
criterion. Since (Hi) contains {iv) as a special case, there remain two implications to 
show. 

{i) ^ (ii) From Sylvester's criterion we know that [£]s,s > for all S* C {1, . . . , n}, < 
5*1 < n. We will show that having a proper principal minor equal to zero will imply 
that zero is not a simple eigenvalue of L under the given assumptions. By eigenvalue 
interlacing [TT] we know that for each symmetric matrix B e R"^" and each S C 
{l,...,n},0<|S'|<n the following inequality holds: 

(B) < Xn.in{Bs,s) 

where A,nin denotes the minimal eigenvalue. Now assume that [i]s,s — for some 
5 C {l,...,n},0 < l^l < n. Then there is i e {1, . . . ,n}\S. Define S = 
Then as seen above = Xmin{Ls,s) > Aminl^g s) — ^^i^ thus [L]g § = and there 
is an eigenvector v oi Lg g for the eigenvalue zero. We will now show how to get an 
eigenvector of L for the eigenvalue zero independent of 1 thus giving a contradiction 
to zero being a simple eigenvalue of L. 



Assume wlog that i = n. Then, because of the zero row sum assumption, L can be 
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decomposed as 



But then 







( 






\ 








L 




Ll 








V 1 




l^Ll ) 




/ 










f 


\ 




L 


-LI 








Lv + 








ii: 








V 




l^Ll ) 






[ - 





(v) (i) If L is not positive semi-definite or has rank smaller than n — 1, there is an 
eigenvalue A < of L and a corresponding eigenvector v = (wi, . . . , w„)-^ independent 
of 1. Define u — v — ti„l and u = {ui, . . . , Un-i) ^ 0. Since u„ = we have 

M^-^-{l,...,n-l},{l,...,n-l}'" = u^Lu 



= v'^Lv - 2vnl^Lv + vll^Ll 
= v'^Lv = Xv^v < 0. 

But this contradicts ..,„_!} being positive definite. 



□ 



3 The Principal Minor Matrix Tree Theorem 

We are looking for properties of a symmetric, zero row sum n x n matrix A C R"^" 
that are necessary for negative semidefiniteness. These properties are best described 
in terms of the structure of A, more precisely, the position and weight of the nonzero 
entries of A. Graph theory provides suitable concepts such as subgraphs, trees and 
cycles to describe this structure. In the light of remark [2] A is negative semidefinite 
iff C{G) — —A is positive semidefinite where G = G{A) is the Coates graph of A. In 
this section, we therefore seek a combinatorial interpretation of the minors of C{G) in 
terms of the subgraphs of G. For S'C{l,...,n},5^0, we define 

J-s ■— {K C E{G) I \K\ — 15*1, each connected component in K contains at least one 
vertex not in S}. 



The definition of J-s implies two features which are characterized by the next lemma. 

Lemma 6. Let K e Ts- Then the following holds: 
(a) K is a forest. 

(h) each tree in K contains exactly one vertex not in S. 
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b 



■ > o 



Figure 1: Examples for J\g. Shown are the sets Ts of a graph G (upper panel) for 
S = {1, 3, 5} (a), S = {1, 2, 3, 4, 5} respectively (b). Panel (b) illustrates the 
relation between Theorem [7] and Kirchhoff's Theorem: For all G and all sets 
S with |5'| = n — 1, J-s is the sum over all spanning trees of G. 



Proof, (a) Assume K is not a forest, i.e. it contains a cycle. Each connected component 
C of K has at most \E{C)\ vertices in S and the component C" containing the cycle 
contains at most |i?(C")| — 1 vertices in S. This gives 

\s\< inc)n5|<(5]|ii;(G)|)-i = |/^|-i = |5|-i 

C is conn. comp. oi K C 

which is a contradiction. 

(b) Assume there is a connected component C of K with more then one vertex outside 
S. Then C" contains at most |i?(C") | — 1 vertices in S and we get the same contradiction 
as in (a). □ 

Theorem 7 (Matrix Tree Theorem for principal minors and undirected graphs). Let 

G he an undirected weighted graph, S C V{G). Then 



[C{G)]s^s = E ^(^) 



where uj{K) = uj{e) 



eeK 



This theorem can be proved using combinatorial arguments as in [5, , but the well 
known proof of the classical matrix tree theorem using the Cauchy-Binet theorem can 
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easily be adapted to this generalization. To keep the paper self-contained we include 
the proof. Its structure is as follows. The Laplacian of the graph is expressed as the 
product of its incidence matrix, its transpose and a weight matrix (Lemma (H]). By the 
Cauchy-Binet theorem the minors of this product are expressed as a sum over products 
of incidence matrix minors (Lemma |5]). Finally these minors are interpreted in terms 
of J^s (Lemma [TO]). 

Lemma 8 (Cauchy-Binet-Theorem for minors |18|). For an nxm matrix D , an mxn 
matrix E and index sets I, J ^ {0, . . . ,n} with |/| = | J| < to we have 

K,J 

KC{l,...,m},\K\ = \I\ 

For the proof of Theorem [7] we first have to introduce some notation: First we impose 
an arbitrary orientation on G. Note that the proofs ofl^andlTOlare independent of this 
orientation since the entries Mie of M occur only in expressions of the form Mi^Mje 
or Mif, + Mjf, for some e £ {1, . . . , \E\}. 

The incidence-matrix M £ {-1,0, l}"^!^^^)! of G is then defined by 

{1 if e starts at i 
— 1 if e ends at i . 
if i ^ e 

Additionally, we define the weight matrix W G as the diagonal matrix 

with Wee — "^(e)- 

The importance of these incidence matrices for our application arises from the following 
two lemmas. 

Lemma 9 (Representation of C{G) by the incidence and weight matrix). 

MVFM^ = C{G) 

Proof. For the off-diagonal entries [MW M'^)ij with i ^ j we have 

(M1¥M^)„ = M,eWeeM,e 
eeEiG)" ' 

Jo ife^{i,j} 
[-a;(e) iie^{i,j} 

which equals C{G)ij. The same holds for the diagonal entries 

{MWM^),,= J2 MlWee= ^ie) = ^G)u. □ 

eGB(G) eeE(G) 
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Figure 2: Illustration of the induction step used in the proof of Lemma 1101 Consider 
the graph G from Fig. [TJ As shown in Fig. db), every tree in Ts^ S = 
{1, 2, 3, 4, 5}, has at least one leaf i in S* and an edge e S £ incident with i. 
Comparison with the sets Tg^^^iy, J-s\{2}j displayed in (a), (b) respectively, 
reveals that for all trees and all i, K\ {e} G Ts\{i}- 



Lemma 10 (Characterization of the incidence matrix minors via .Fs)- 

\[M]s.k\ = r If^l for allSCV,KC E, \K\ = \S\ 

ifK^Ts 



Proof. {K ^ Ts ^ [M]sj< ^ 0) li K ^ Ts, there is a connected component C of 
that has no vertex outside 5. Then the rows of Ms.k corresponding to these vertices, 
sum up to 0. 

{K e J-s ^ [M]s,K — 1) We prove this direction by induction on \S\. If |5| — 1 
and K £ J-s then K is a single edge between a vertex in S and one outside S, so 
[-^^]s',i<'| = 1- If 15*1 > 1 then by Lemma [S] there is a tree in K, which has at least one 
leaf i in 5 and an edge e € K incident with i which implies that K\{e} G J's\{i}- Since 
the row of Ms.k corresponding to i contains only one nonzero element, by expansion of 
the determinant along this row |[M]5^if | — \ [M]g\^iy ]^\^^^y \ which implies the assertion 
by induction. □ 



We are now ready to proof Theorem [71 
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Proof of Theorem^ Using Lemma 9, Lemma 8 and Lemma 10, we get 
[CiG)]s,s^[MWM^]s,s^ J2 [MW]s,k[M^]k,s 



KCE 
\K\ = \S\ 



emi,^ n ^(^) 

K eeK 

E ^w- ° 



4 Cutting a Forest 

We now have a convenient characterization of the minors of a zero-row-sum symmetric 
matrix L G M"^" in terms of forests contained in G = Q{—L). To show that a given 
structure of G is incompatible with the positive semi-definiteness of L, one needs 
to show that the inequalities [i]s,s > 0, S ^ {l,...,n}, cannot simultaneously be 
fulfilled. 

However, if one can show that e.g. X^sWs.s — one may dispense with the explicit 
calculation of any of the inequalities as well as with the compatibility verification. In 
this section, we prove a combinatorial identity that is central for this type of calcula- 
tion. 

Let Vi, V2 be a partition of the vertex set of G into two non-empty parts. For -B C Vi 
define 

Es = {D C £: n E{Vi,V2) I D is a forest and 

contains exactly one edge e for each i £ B with i e e}, 
Tb = {D CEr\E{Vi,Vi) \ \D\ = \Vi\ - \BlD is a forest and 

each tree contains exactly one vertex of i?}. 

Lemma 11 (Cutting a forest). Let C CVi. Then 

^V^\C^ U {AiJA'\A£EB,A'£TcuB}- 
BCVi\C 

Proof. (C). For K £ J-Vi\c we clearly have K f] E{Vi,V2) £ ^b where B is the set of 
all indices of Vi incident to an edge of K between Vi and ¥2- By definition of JVi\c 
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we also have K n E{V2,V2) = and therefore \K n EiVi, Vi)\ = \Vi \C\~ \B\. This 

gives if ns(Fi,yi)ercuB. 

(3). On the other hand let Ks G'Sb and Kt £ Tsuc- Then each tree in K := i^sUi^T 
contains exactly one vertex not in 14 \ C and \K\ = - |C U S| + \B\ = \Vi \ C\. 
Therefore K e Tvi\c- □ 

This decomposition of [C{G)]y^\c .Vi\c ahows us to prove the following remarkable 
identity. 



Theorem 12. 



^ (-l)l^l..(Sc)[/:(G)]v,\c,yAc = 



Proof. From Lemma [TT] we can derive 



J2 E ^(^5)^(/^t) 



i^T ^ T'b U C 



E f E -(^-)) ( E 



uj{Kt) 

BCVi\C XKsESb / KKtGTbuc / 

= E w(SB)a;(rBuc), 

BCVi\C 

E (-l)""^(Sc)[AG)]v^,\c,yAc- E (-l)'^'^(^^) E ^(SB)c.(rBuc) 
ccvi ccvi Bcyi\c 

= E (-l)''^'^(Ssuc)'^(7suc) 

CCVi 
BGVi\C 



ACVi 
CCA 



E ^(Sa)c.(T^)E(-1)""=0, 



ACVi CCA 



where we use the fact that 



CCA k=l ^ ^ 
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5 Positive Spanning Trees 



A weighted graph is said to be positive or negative, if each of its edge weights is positive 
or negative, resp. 

Lemma 13. Let G be a weighted connected graph. Then G has a positive spanning 
tree iff it has no negative cut. 

Proof. Take an arbitrary cut of G given by the partition Vi, V2 of If G has a positive 
spanning tree, then there must be an edge from Vi to V2 with positive weight, so the 
cut is not negative. If on the other hand there is no negative cut, remove all negative 
edges from G. Then each cut still contains an edge and therefore the remaining graph 
contains a spanning tree with all edge weights positive. □ 

From now on let A be a real symmetric negative semi-definite matrix with zero row 
sum. 

Theorem 14. If G = G{A) is connected then G{A) has a positive spanning tree. 

Proof. Assume there is no positive spanning tree. By the above theorem this implies 
the existence of a negative cut given by the partition Vi, V2 of V{G{A)), that is Ve 6 
VixV2nE: w(e) < 0. For i € V{g{A)) define z, = J2j(,v, and z'^ = Y.jev, • 
Since A has zero row sum, Zi + z,[ = und by the negativity of the above cut Vi G Vi : 
z[ < 0, 3i e Vi : z ■ < 0. Finally consider w e M" with 




1 if i e 
if i e V^2 ■ 



Then v'^Av = J2i jeVixVi ~ SieVi -^^ ^ ^' thus A is not negative semi-definite, a 
contradiction. □ 

Remark 15. For matrices of rank n—\, one can also derive this result from Theorem 
Since Vi,V2 define a negative cut, (— l)l'^la;(I]c') non-negative for all C and 
positive iff all vertices in C are incident to an edge between Vi and V2 . Since according 
to Theorem\^all minors C{G)s.s with \S\ < n are positive if our system is stable, the 
left-hand-side of the identity in Theorem is positive, whereas the right-hand- side is 
zero. 

Remark 16. IfQ{A) has disconnected components, every such component has a pos- 
itive spanning tree. 

We now apply Theorem [T3] and Theorem [7] to derive an upper bound on the absolut 
weight on negative edges in induced lines. 
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Theorem 17. Let H be an induced line ofG{A) = {V,E,uj) with at least two edges, 
i.e. a connected induced subgraph for which all but two vertices have degree two in G{A). 
Then H contains at most one negative edge e and if it contains one, the absolute value 
of its weight is bounded from above by 

IlaeE{H)\{e} ^(a) _ 1 

J2beE{H)\{e} nQe_E(_fi-)\{e,b} '^('^) Yl^aeE{H)\{e} ZJ^^ 

Proof Let T be a spanning tree of GiA). Then \E{H) \ E{T)\ < 1, hence H contains 
at most one negative edge by Theorem [T31 By definition of an induced hne, H is a 
tree with two leaves. Let {zi, . . . , i\H\} — ^(H) be an enumeration of the vertices of 
H such that E{H) = {{ik,ik+i) | fc S {1, ■ • • , \H\ — 1}}. The bound can immediately 
be derived from Theorem [7] applied to 5 = {«2, . . . , Notice that the elements 

of !Fs are precisely the subsets of E{H) of size \H\ — 2. The bound on the weight 
of a negative edge is then derived by using the non-negativity of the principal minor 
[A]s.,s- □ 

Every connected subgraph of an induced line is again an induced line. Hence for each 
induced line in which a negative edge is contained, Theorem [T7] gives a bound on the 
weight of this edge. If one such line is contained in another, the longer one yields a 
sharper bound. 

Example 18. We now apply the results of this section to the Kuramoto model defined 
in Eq. ([1]) and reveal common properties of the configuration of all possible phase- 
locked systems. Consider Theorem^T^with the edge weights graph G = G(A) given by 
Aij = Bij cos{(fji). Since all coupling constants Bij > 0, .stability of the phase-locked 
state requires that each component of the physical coupling network has a spanning tree 
of oscillators obeying \ipji \ < ^. 

Theorem \ 1 7| implies that .stability of the phase-locked state can only occur if along each 
induced line H of G there are at most two neighbouring oscillators with absolute phase- 
difference greater or equal ^. If such a pair (i,j) of oscillators exists, the phase-locked 
state can only be stable if the absolut value of the reciprocal of cos{ipji) is larger then the 
coupling constant Bij times the sum of the reciprocals of all other coupling constants 
Bki in H. 



6 Conclusion 

In the present paper we used an extension of Sylvester's criterion to formulate meso- 
scale obstructions to stability of one-dimensional center manifolds, in particular, of 
the center manifold of phase-locked solutions for the Kuramoto model on complex 
networks. We obtained a number of necessary conditions for stability which restrict 
the position and weight of edges of the system's Coates graph. Among other results. 
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this investigation provides a rigorous proof for the so-called positive-spanning-tree 
criterion ^ that states that stability in the Kuramoto model requires the existence of 
a spanning tree of edges associated with positive elemets of the system's Jacobian. 

In network dynamics, system-level implications of meso-scale structures are just now 
intensively studied [9l [1] . While efforts were long concentrated on the implications of 
local topological properties (such as the degree distribution), and global topological 
properties (such as the clustering coefficient), recent results [H [151 E] indicate that 
also meso-scale properties can crucially affect synchronization. We believe that the 
results obtained here will be valuable for network dynamics, because they allow to 
extract meso-scale structural properties for the dynamics of networks. For a minor of 
given size fc, already an application of Sylvester's criterion reveals necessary conditions 
for stability that restrict the structure of motifs containing k vertices. The conditions 
obtained in this way are only necessary and not sufficient because instability could 
still occur due to instabilities in larger or smaller motifs. More complex necessary 
conditions for stability can be formulated on the meso-scale. E.g. for induced lines 
(Theorem [T7| , an obstruction to stability are lines with more then one negative edge 
or those with one negative edge whose weight is below a certain threshold. 

In the future this approach may allow to attribute dynamical characteristics such as 
the onset of oscillations or the loss of synchronization to a specific meso-scale structure 
in the system. Such insights are presently actively sought in network dynamics |13j . 
Further, they may hold the key to important real world applications such as the 
operation of smart grids, future adaptive power grids, which will have to maintain 
synchronization in very different demand situations. 
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